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ABSTRACT 

The formation and evolution of low-mass stars within dense halos of dark matter (DM) leads to 
evolution scenarios quite different from the classical stellar evolution. As a result of our detailed 
numerical work, we describe these new scenarios for a range of DM densities on the host halo, a range 
of scattering cross sections of the DM particles considered, and for stellar masses from 0.7 to 3 M©. 
For the first time, we also computed the evolution of young low- mass stars in their Hayashi track in the 
pre-main sequence phase and found that, for high DM densities, these stars stop their gravitational 
collapse before reaching the main sequence, in agreement with similar studies on first stars. Such stars 
remain indefinitely in an equilibrium state with lower effective temperatures (|ATe//| > 10^ K for a 
star of one solar mass), the annihilation of captured DM particles in their core being the only source 
of energy. In the case of lower DM densities, these proto-stars continue their collapse and progress 
through the main sequence burning hydrogen at a lower rate. A star of 1 Mq will spend a time greater 
than the current age of the universe consuming all the hydrogen in its core if it evolves in a halo with 
DM density = 10^ GeV cm~^. We also show the strong dependence of the effective temperature 
and luminosity of these stars on the characteristics of the DM particles and how this can be used as 
an alternative method for DM research. 

Subject headings: dark matter - Galaxy: center - Hertzsprung-Russell diagram - stars: evolution - 
stars: formation - stars: interiors 



1. INTRODUCTION 

Modern observational Cosmology has revealed a more 
complex and unknown Universe than previously ex- 
pected. The current cosmological observations of the 
baryon acoustic oscillations, distance measurements by 
means of Type la super novae, cosmic microwave back- 
ground and primordial light elements abundances, are 
all observations that confirm that the standard cosmo- 
logical model is a good description of our Universe un- 
dergoing an accelerated expansion fSpergel et al.l [2007[ 
iKomatsu et al. 2009). We now know that only about 5 
% of all matter in the Universe is regular visible matter, 
usually known as baryons. The rest is very likely to be 
in the form of unknown particles that have not yet been 
detected, even though we are able to observe its effect on 
the formation of cosmological structures, usually known 
to us as dark matter (DM) and dark energy. The Uni- 
verse is composed of 20% of DM, particles that do not 
interact with the electromagnetic field, but whose pres- 
ence can be inferred from gravitational effects on visi- 
ble matter. The remaining 75% is known as dark en- 
ergy. This component is very homogeneous and not very 
dense, and is known to have interacted through any of 
the fundamental forces other than gravity. In this paper 
we investigate the impact of the DM component on the 
evolution of stars. 

There is an abundant number of particle physics mod- 
els providing candidates for DM particles. Among oth- 
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ers, the supersymmetric standard model (SUSY) is one of 
the best studied candidates for physics beyond the Stan- 
dard Model (Bertone et al. 2005). Dark matter particles 
must be massive particles, electrically neutral and, very 
likely, non-coloured. A lightest supersymmetric parti- 
cle (LSP), such as the neutralino, is the favourite SUSY 
particle among the particle physics community. The LSP 
belongs to a generic family of neutral massive particles, 
usually known as Weakly Interacting Massive Particles 
(WIMPs), that are by definition the best DM particle 
candidates. WIMPs interact through gravity and pos- 
sibly through the weak nuclear force, and very likely 
through no other interactions stronger than the weak 
force. 

Despite their weak interactions, WIMPs can lead to 
significant changes in the formation and evolution of 
stars, provided they have a sizeable scattering cross sec- 
tion with baryons. Concerning WIMPs' interaction with 
baryons, the consequence is twofold: WIMPs can affect 
a star by annihilating among themselves into standard 
model particles in its core, providing a source of en- 
ergy additional to standard nuclear energy (jSalati fc Silkl 
119891 ). A second way WIMPS can influence stellar struc- 
ture is by providing an additional mechanism of heat 
transport inside the star (Bouquet & Salati 1989a). This 
can reduce the local temperature gradient, potentially 
inhibiting convection an d enhancing the pul sation of 
horizontal branch stars (|Dearborn et al.l Il990f ). In the 
case of the Sun, the current values favoured by WIMP 
self-annihilation cross section, WIMP-nucleus scattering 
cross sections and local DM density indicate that this 
effect is not significant for the evolution of our star. 

The study of the effects of DM particles accretion on 
stars is an alternative approach to investigate the prop- 
erties of such particles. These effects were first studied 
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in the Sun by ISpergel fc PressI (IT985I) and later on by 
Lopes et al.l (l2QQ2af ) and lBottino et alJ ()2QQ2l). Recently, 
Moskalenko fc Wail ^200?) for white dwarfs. [Scott et al.l 
(|2QQ7f ) for low-mass stars, and iSpolvar et al.l (j2QQ8f ) for 
the first generation of stars, showed that, if embedded 
in halos with high DM densities, stars can be fuelled 
only by the energy from DM annihilation. Many au- 
thors confirmed and developed these results. DM capture 
and annihilat ion was studied for white d warfs and neu- 
tron s tars by Bertone fc FairbairnI (|2Q Q8^ and Kouvaris" 
(l2008h. Th e' same wa s done for fi r st star s by locco (2008), 
iTaoso et al.l (|2008f ). IFreese et al.l (j2008af ). and lYoon et al.l 
()2008f ). and their implications for cosmic reionization 
and the pair pro duction su pernov a rate were studied by 
iSchleicher et al.l (j2009f ) and lloccol (j2QQ9f ). respectively. In 
the case of first stars formed at the center of DM mini- 
halos, DM annihilation may become the primary source 
of energy during the adiabatic contraction regime of the 
DM halo, counteracting the gravita tional collapse even - 
before DM ca pture becomes efficient ("Spolvar e t al.|[2QQ8l: 
Freese et al.ir2 QQ8b: Freese et al. 2008c; Natar a^an et al.l 
20091 : iRip amonti et al.ii200 9). When the star runs out of 
its or iginal DM, which takes about a million years (Freese 
|2009), this stalling phase ends and the collapse contin- 
ues until cap ture replenishes the fuel for DM annihilation 
(jlocco et all [2008). 

Regarding low-mass stars, iScott et al.l ()2009f ) carried 
out an extensive study on how the stars located at the 
galactic center may be affected by the dense DM den- 
sities which are expected to exist there. They assumed 
that these stars were already formed in a scenario with- 
out DM. Therefore, their evolution on a DM halo was 
considered from the ZAMS. In addition to that, we also 
took into account the capture of DM p articles from the 
collapse of the protostar, as l occo et al. (2008) did for 
first stars. The influence of DM in that early stage may 
have dramatic consequences in the forming star, as will 
be shown. 

The paper is organised as follows: in section [2] there 
is a short description of the stellar evolution code used 
here. The process of accretion of annihilating DM par- 
ticles in stars is briefly reviewed in section [3l The dif- 
ferent evolution paths found for low-mass stars evolving 
within DM halos are discussed in section HI In section [5] 
we discuss the implication of our results to observational 
stellar astrophysics. Finally, in section [6l we present a 
brief summary of the results and discuss the implication 
of such results for DM research. 

2. THE STELLAR EVOLUTION CODE 

The models computed in this work were made using 
the stehar evolution code CESAM (Morel 1997). The 
CESAM code is a consistent set of programs and rou- 
tines which performs calculations of ID quasi-static stel- 
lar evolution including diffusion, rotation and mass loss. 
CESAM computes structure equilibrium models by col- 
location method based on piecewise polynomials approx- 
imations projected on their B-spline basis, the evolu- 
tion of the chemical composition is solved by stiffly sta- 
ble schemes of orders up to four, the solution of the 
diffusion equ ation employs the Petrov-Galerkin scheme 
(jMorell [l997[ ). The code determines the evolution of a 
star by the integration of the set of conservation struc- 
ture equations, coupled with a set of nuclear reactions 
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describing the nucleosynthesis of chemical elements and 
the production of energy, using an adaptive space-time 
mesh, i.e., the code chooses an optimal step in time and 
space by computing the rate of variation of the equi- 
librium quantities. CESAM allows calculations of stellar 
models with various physical assumptions, physical data, 
external boundary conditions, numerical methods and 
numerical accuracy; in this work the accuracy is set to 
10 This code has the ability to compute the evolution 
of stars from the pre-main sequence up until the begin- 
ning of the ^He burning cycle for various stellar masses, 
and for a range of metallicities (0.0004< Z <0.04, in 
the present configuration) . The initial chemic al abun- 
dances are the solar ones (|Asplund et al.|[2QQ5l ) and the 
initial metallicity, unless stated otherwise, is assumed to 
be Z=0.014, similar to that of the Sun. 

The microscopic physics, including equation of state, 
opacities, nucleosynthesis, microscopic diffusion and 
chemical abundances, is very refined in this stellar evo- 
lution code, since this part of the code was tested in the 
case of the Sun against helioseismic data. The CESAM 
evolution code is well established in the solar and stel- 
lar phys ics comunit y, being u sed either to comp ute solar 
models (jCouvidat et al.ll2QQ3l : iBerthomieu et al.|[l993t ) or 
models of other stars (anions o thers: PMS(5Scuti star 
V346 Ori ([Bernabei et all I2009D . /3Chepei star u Eri- 
dani jS uarez et al. 2009), 0.8-8M,7, PMS and MS stars 
(jMargues et al. 2008), Vega-like stars aPsA, /3Leo, /3Pic, 
eEii and rCet (iDi Folco et al.i 2004), aCMi star Pro- 
cyon A (Kervella et al.ll2004f). solar-like 0.8 -1.4Mr7^ stars 
and P MS stars dPiau fc Turck-Chiezell2002l ). aCen AfcB 
stars (|Thevenin et al.ll2002D and MS stars with convec- 
tive overshooting (|Audard et al.l [19951 )). 

We modified this code to take into account the im- 
pact of DM particles in the evolution of low-mass stars. 
In particular, we have already used this code to study 
the evolution of th e Sun within a halo o f dark matter 
(jLopes et al.l 1200 2 al ibi: iLopes fc Siffll2002D : we predicted 
the results of two groups of observables: solar neutri- 
nos and helioseismology data (including the sound speed 
profile). 



3. ACCUMULATION OF DARK MATTER PARTICLES 
INSIDE STARS 

WIMPs travel through stars, where they experience 
scattering off the nuclei that they encounter in the stel- 
lar cores. Although most WIMPs travel right through 
the star without suffering any type of interaction, some 
of them will scatter off nuclei losing energy. If they lose 
enough energy, they won't be able to escape the gravita- 
tional field of the star anymore. The number of WIMPs 
trapped inside the star is measured by the capture rate 
C^. The capture rate is proportional to the WIMP scat- 
tering cross section off nuclei and the dark matter 
density of the halo p^. It is inversely proportional to 
the WIMP mass and to the WIMP dispersion veloc- 
ity v^. Capture r ates o f WIMPs were first calcul ated by 
iPress fc Spergell (|1985f ) in the case of the Sun, by iGouldl 
(|1987l ) for generi c massive bodies and for the earth in 
particular, and by "Bouquet fc Salati' Cl989al) for main se- 
quence stars. We computed the capture rate using equa- 
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tion 2.31 of lGouldl (UqSTI) : 
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where Vesc = ^JlGM^j is the escape velocity, and 
are the mass and radius of the star, is the proton 
mass, and xi and Ai are, respectively, the mass fraction 
and the number of nucleons of element i. We included the 
contributions from 13 nuclei: H, "^He, ^^C, ^^n, i^q^ 2jj^ 
^He, ^Li, ^Be, ^^C, ^^N, ^^O and ^Be, the abundances 
of which are followed by our code. For all of them, we 
computed their spin-independent (SI) interaction with 
WIMPs (coherent scattering), which scales as the fourth 
power of the nucleus mass number, leading to a scatter- 
ing cross section Oi = Afa^^si- For hydrogen, we also 
took into account the contribution of the spin-dependent 
(SD) interaction (incoherent scattering), which lead to 
an = Ajja^^si (^x^sd- The last term in equation [TJ 
is the product of the suppression factors due to the mo- 
tion of the star through the halo ^^(oo), the mismatch 
between WIMP and nuclei masses, and the dimensionless 
gravitational potential averaged over the star < (j) > , all 
present in the original expression of lG^ ([1983). We 
evaluated these factors for a WIMP mass m^ = 100 GeV, 
a velocity of the star v^ = 220 km s~^ and a WIMP 
dispersion velocity = 270 km s~^, and found that, 
for these values and for all stars studied here, we can 
approxin iate C to the order of un i ty, similarly to other 
authors ("Moskalenko & WaJ 120071 : iBertone fc FairbairnI 
120081 : iFree se et al. 2008a). This approximation retains 
the main aspects of WIMPs' capture while it simplifies 
the calculations and reduces the time of computation. 
Furthermore, the difference in values between this cap- 
ture rate and others (see table [1]) is relatively small. We 
confirmed that the overall conclusions of this paper will 
not be affected by the accuracy of this capture expres- 
sion. 

The values of the scattering cross sections used 
our simulations, ct^^sd from 10~^^ cm^ to 10~^^ cm 
and cr^^si = 10~^^ cm^, are consistent with the 
experimental bounds given by d i rect detection ex- 
periments ( C DMS (jAkerib et al.l JSoog), XENONIO 
(jAnde et al ] |2QQ8bD, NAIAD ( |A 
PICASSO (iBarnabe-Heider et al.l 
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(Behnke et alT 
Ahmed et al. 
CRESST-II 



ner et al.l [2005), 
20051) . COUPP 



for a^s D. and CDMS-I I 
2009), XENONIO (lAnde et al.l [2008a), 
(Anglohe?llS] [20051) and EDELWEISS 
([Sanglard et al.l 12005) for a^^si)- Only one indirect 
detection experiment (Super-Kamiokande (|Desai et al.l 
[2004)) predicted an upper limit for (Jx.sd below 10~^^ 
cm^. The bounds on the WIMP-nucleon scattering cross 
sections are much less constraining for SD interactions 
than for SI, due to the presence of the A^ factor in the 
latter type of interactions. For this reason, the WIMP 
capture rate was always dominated by SD scattering in 
all our computations. 

After being captured, WIMPs will then sink to the core 
of the star, where they can annihilate with one another at 
an annihilation rate F^. The total number of WIMPs A^^ 
inside the star is then determined by the balance between 
the capture rate and the annihilation rate. Therefore, the 
evolution of the number of WIMPs inside the star over 
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TABLE 1 

Comparison of different capture rates in the literature, 

considering different stellar masses and dm halo 
DENSITIES p^. In all CASES: M^ = 100 GeV, (7x,sd = l^~^^cy]? . 

=0.01, 



The metallicity of the stars goes as follows: 
*:Z=0 AND *:Z=0.018. 



t:Z= 



time is 



(2) 



The capture and annihilation of WIMPs in the core of 
a star is a very efficient process; capture and annihi- 
lation processes_balaji^ each other out in around 100 
years (|Salati fc Silklll989l ) and the system rapidly comes 
to equilibrium, = or = 2r^. The time scale 
for the steady state to be reached will be much inferior 
than the typical evolutionary time scale of a main se- 
quence star. WIMPs will get distributed in the core of 
the star very rapidly following an approximately ther- 
mal int ernal distributio n with a characteristic radius 
= 3i<iTc/ 27TGpcm^ , where Tc and pc are the cen- 
tral temperature and density of regular baryonic matter 
inside the star, and G and hz are the Newton and Boltz- 
man constants. The density number distribution is 

given by n^{r) = A/"^7r~^ ^^r~^ (i Dearborn et al.l 

1990, l£l9llKaDlan et al.llT99il ). 

Dark matter particles can transport energy by scatter- 
ing off nuclei inside the star, thus con stituting an alter- 
native mechanism of energy transport CBouQuet fc Salatil 
.1989b.) . This energy transport mechanism by WIMPs 
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is implemented in o ur stellar evolution code (for a de- 
tailed descripton, see lLopes et alJ ()2QQ2bl )). although the 
contribution of such process to the evolution of the star 
is negligible compared with other effects. More signifi- 
cantly, DM particles in the stellar core provide an extra 
source of energy (Salatijfe Silk 1989). The energy gener- 
ation rate due to pair annihilation of DM particles is 
given by 



^x(0 = fx ri\{r)p{ry 



< CFaV >, 



(3) 



in units of energy per mass per time (in cgs: erg g~^ s~^). 
It follows that every pair of DM particles captured into 
the star is instantly converted into additional luminos- 
ity. We assume that all products of DM annihilation, 
except neutrinos, interact either by electromagnetic or 
nuclear strong forces with the core nuclei, so they have 
short mean free paths inside the star and these parti- 
cles rapidly reach the thermal equilibrium. We chose 
the coefficient = | x 2 to take into account that | 
of the energy is lost in the form of neutrinos that es- 
cape from the star, and that each annihilation involves 
2 DM particles, due to the assumption that the neu- 
tralino is a Majorana particle, i.e. it is its own anti- 
particle. The coefficient has different values in the 
literature. Its former factor, which quantifies the en- 
ergy that remains inside the star, could be underesti- 
mated in o u r wor k: following the recent simulations of 
iScott et all (j2QQ9[ ). the energy loss could be as low as 
10% of the total energy from DM annihilations. Our 
choice, more conservative, is in agreenient w i th other au- 
thors (jFreese et al.ll2QQ8cl : llocco et al.ll2QQ8l : lYoon et al.l 
[2008h . In our simulations we assumed the annihilation 
cross section to be < (JaV >= 3 • 10~^^ cm^ s~^, a value 
that is fixed by the relic density through the following 
appr oximation: Vt^h'^ ^ 3 • 1Q ~ ^^ cm^ s~-^/ < crgV > 
(e.g., iScherrer fc Turnerl (jl986f ): iSrednicki etaU (jigSSD : 
iBertone et all (j2QQ5f )). 

4. EVOLUTION OF LOW-MASS STARS WITHIN A DARK 
MATTER HALO 

We have implemented the effects of the annihilation of 
DM particles in our stellar evolution code and followed 
the evolution of low-mass stars since their proto-star 
phase and throughout the main sequence phase. These 
stars may experience dramatic changes on their evolution 
depending upon the amount of DM the star accumulates 
in its interior. The accretion of DM depends mainly on 
the particle halo density p^, and also on the WIMP- 
nucleus spin-dependent scattering cross section cr^^sD- 
The more accretion of DM particles happens inside the 
core of the star, the more energy is produced by WIMP 
pair annihilation. The existence of this new source of 
energy leads to significantly different scenarios of stellar 
evolution. Figure [T] shows the contribution of the differ- 
ent energy sources to the total energy generation rate, 
^T, of a star of 1 Mq. The evolution of the star depends 
on the balance between DM energy rate, e^, the ther- 
monuclear energy rate produced by the pp chain, e^p, the 
thermonuclear energy rate produced by the CNO cycle, 
^CNO^ and the gravitational energy rate, Sgravi produced 
by the gravitational contraction of the star. Depending 
upon the amount of DM present in the host halo, we 
found that stars can experiment quite different evolution 
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Fig. 1. — Energy rates during the evolution of a 1 Mq star on the 
Classical, Weak and Strong cases. The densities of dark matter in 
the halo are, respectively: 0, 3- 10^ and 3- IQ-'^-'^ GeV cm~'^. Curves 
are as follows: et (light grey dashed line), (red continuous line), 
Egrav (bluc dashcd line), Spp (green long dashed line), ecNO (pink 
dotted line). The DM tialo is assumed to be composed by WIMPs 
with mass = 100 GeV, spin-dependent scattering cross section, 
crx,SD = 10" 
3 • 10-26 cm3 



cm^, and annihilation cross section < aaV >- 



paths, which we classified in three distinct cases: Weak, 
Intermediate and Strong case scenarios. 

4.1. Scenarios of stellar evolution within DM 

4.1.1. Weak case scenario: Slowly evolving stars 

Normal stars are self-gravitating systems that most of 
the time are experimenting a gravitational contraction, 
leading to an increase in the temperature inside their 
cores. The gravitational collapse is stopped by an ad- 
ditional source of energy, such as thermonuclear energy 
produced by the pp chain or CNO cycle in stars on the 
main sequence phase. Nevertheless, stars evolving in DM 
halos can experiment a quite different scenario of evolu- 
tion. For stars evolving within halos with low DM den- 
sity p^, the energy from WIMPs' annihilation is a com- 
plementary source to the thermonuclear energy (see Fig- 
ure [TJb), that slightly delays the gravitational collapse, 
slowing down the arrival of the Hydrogen burning phase. 
The equilibrium is reached at a lower central tempera- 
ture than that of the classical evolution case, leading to a 
smaller rate of energy produced by thermonuclear reac- 
tions £pp + sc NO (stars will evolve in the weak scenario if 
their thermonuclear energy accounts for more than 10% 
of the total energy in the beginning of the MS). There- 
fore, the time that a star spends in the main sequence 
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Fig. 2. — Time spent in main sequence (from Egrav < 1%8t to 
Xc < 0.001) by stars of masses 1, 1.5, 2 and 3 M© embedded in 
dark matter halos of different densities. The DM particles are as 
described in Figure [T] 

phase is enlarged respect to the classical evolution sce- 
nario (see Figure [2]). The more massive the star, the 
more DM will be necessary to produce the same effects. 
A star of IM© will stay in the MS for a time greater than 
the current age of the universe if it evolves in a DM halo 
of = 10^ GeV cm~^ and cr^^sD = 10~^^ cm^, while 
a star of 2 Mq evolving in the same halo wil l not be af- 
fected. In the case of one solar mass stars, iScott et al.l 
(12009) obtained the same extension in the main-sequence 
lifetime for almost identical DM densities on the host 
halo. On the other hand, for greater masses our results 
are more conservative, due to the lower WIMP capture 
rates obtained for > 1 Mq. This evolution scenario 
is qu alitatively s imila r to that predicted for Pop III stars 
bv lTaoso et al.l (j2008f ). 

To grasp the role of the metallicity, we computed 
models with metallicities from Z=0.0004 to Z=0.04 and 
found that the main differences in the stellar evolu- 
tion are those already expected in the classical picture; 
stars with higher metallicities have lower thermonuclear 
energy production rates and, therefore, extended main 
seque nce lifetimes (Schaller e t al.l 119921 : iSchaerer et al.l 
119931 : iLeieune fc Sc haerer 2001. Table [2] shows the 
energy rates Spp and for stars with metallicities 
Z=0.0004, 0.001, 0.02, and 0.04 that evolve in a halo with 
a DM density = 10^ GeV cm~^. Even though is 
lower in high metallicity stars (due to a lower capture 
because of the smaller hydrogen mass fraction), its per- 
centage over the total energy is higher given the strong 
reduction in Spp. Consequently, stars with higher ini- 
tial metallicities will be affected by the energy from DM 
annihilation at slightly lower DM densities. 

The competition between the nuclear burning and the 
energy from WIMP annihilation leads to another impor- 
tant change respect to the classical scenario. The path 
that these stars follow on the HR diagram may be sig- 
nificantly altered if there is enough DM in their interior. 
This can be seen in Figure [H where we plotted the tracks 
of a star of 1.5 Mq evolving in halos without DM and 
with a DM density p^ = 10^° GeV cm~^. In the latter 
case, the rates of energy production were, in the beggin- 
ing of the MS: £^ 42%, Spp 55% and scno - 3%. 
In this case, the contribution of rapidly starts to com- 
pensate the collapse of the proto-star, stopping it when 
the star has a larger radius than that of the classical sce- 



0.0004 0.7584 

0.001 0.756 

0.02 0.680 

0.04 0.620 



0.2412 
0.243 
0.300 
0.340 



0.8 (12%) 
0.7 (12%) 
0.6 (23%) 
0.6 (25%) 



5.8 (88%) 
5.4 (88%) 
2.1 (77%) 
1.7 (74%) 



TABLE 2 

Energy rates (and its percentage over the total energy) 

FOR STARS OF 1 Mq WITH DIFFERENT INITIAL METALLICITIES, 
EVOLVING IN A HALO WITH A DM DENSITY = 10^ GeV CM~^, 
AT AN AGE SUCH THAT THEIR CENTRAL HYDROGEN MASS FRACTION 
IS Xc = 0.60. 
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Fig. 3. — HR diagram of the formation and evolution through 
the main sequence of a star of 1.5 embedded in halos without 
DM (grey line) and with a DM halo density p^ = 10^° GeV cm"^ 
(red line). The DM particles are as described in Figure [T] 



nario. Consequently, the ZAMS position of these stars 
shifts to lower effective temperatures. 

The internal structure of the star also experiences some 
important changes. In the classical scenario, a Sun-like 
star of one solar mass will develop a small convective core 
within the radiative interior, which will disappear when 
the star reaches the main sequence phase. Then, this star 
will develop a radiative core and a convective envelope 
on the outer layers. Alternatively, in this new scenario, 
the central convective zone remains for a longer period 
during the evolution of the star, because the extra lumi- 
nosity amount produced by the WIMPs annihilation re- 
quires a more efficient mechanism to evacuate this energy 
from the stellar core. Typically, in a halo of DM parti- 
cles with p^ = 10^ GeV cm~^ and cr^^sD = 10~^^ cm^, 
a star of 1 M© develops a convective core with a radius 
of approximately 0.05 R^. This radius decreases with 
time, until it disappears completely at an age of 6 Gyr. 
Similarly to the Sun, this star conserves its convective 
external region throughout its evolution. In the case of 
the present Sun, the external convective region is located 
above 0.7 R^. The thickness of this external convective 
layer grows with p^, and it is also conserved in more 
dense DM halos, as in the case discussed in the next 
section, the intermediate case scenario. 
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4.1.2. Intermediate case scenario: Convective-radiative 
^Jrozen" stars 

As the ambient density of DM particles increases, 
the capture rate inside the star increases too. As a con- 
sequence, the energy source resulting from the annihi- 
lating DM particles eventually starts to compensate the 
gravitational energy source, thus keeping the tempera- 
ture of different core regions below the threshold needed 
to start the thermonuclear reactions. In this scenario, 
WIMP pair annihilation becomes the only source of the 
star's luminosity and is high enough to stop the grav- 
itational collapse at an early stage. Stars in this scenario 
can live forever without the production of thermonuclear 
energy. 

Since WIMP annihilation occurs in a more centralized 
region than the nuclear one, at least for very massive DM 
particles, the temperature gradient is much steeper in the 
core of the star than it would be otherwise, and the star 
has the conditions to maintain a convective core for the 
rest of its life. Outside the convective core, less energy 
is generated per unit volume than if the nuclear burning 
was proceeding normally, so the temperature gradient is 
smaller. The actual temperature is lower than in a nor- 
mal star, however, it remains high enough to prevent any 
major increase in opacity, ensuring that energy transport 
in the region above the core remains radiative. The en- 
ergy from the core is easily transmitted through this ra- 
diative region to the surface of the star. The external 
envelope of this star is very much similar to a typical 
young Sun. In the outer layers, the star develops a con- 
vective region; due to the rapid temperature drop, some 
chemical elements such as oxygen, carbon and nitrogen, 
fully ionised in the interior, are partially ionised in the 
most external layers. This increases significantly the ra- 
diative opacity, and makes convection the only efficient 
mechanism of energy transport towards the surface. 

4.1.3. Strong case scenario: Fully convective ^Jrozen" stars 

At high enough WIMP capture rates, the energy pro- 
duced by the annihilation of DM is sufficient to fully 
compensate the gravitational energy (see Figure [He). 
The star's gravitational collapse stops before reaching 
enough central temperature to begin nuclear fusion, as 
in the intermediate case. This equilibrium is reached 
quite early in the formation of these stars, depending 
upon the value of , as illustrated in Figure HI This 
Figure shows solid grey lines corresponding to the clas- 
sical evolution tracks of stars from 0.7 to 3 M©, along 
with the position on the HR diagram of the early equi- 
librium states reached by these stars, considering forma- 
tion scenarios with different WIMP densities p^. These 
cases are strongly reminiscent of the Hayashi track which 
young stars follow when travelling along the proto-star 
phase onto the final stages of their formation towards 
the main sequence phase. If no extra source of DM en- 
ergy existed, these stars would shrink in size within the 
Kelvin-Helmholtz timescale as they radiate away grav- 
itational energy. Alternatively, a constant energy gen- 
eration in the core by WIMP pair annihilation creates 
stars that can in principle remain in the same position 
in the HR diagram for an arbitrarely long time. A star 
of one solar mass evolving within a halo of DM particles 
with a scattering cross section crx,SD = 10~^^ cm^ and 



J. & Lopes, I. 

3.5 r- 

3 - 



2.5 - 




4.2 4.1 4 3.9 3.8 3.7 3.6 3.5 



log ( Te// / K ) 

Fig. 4. — Stationary states reached by stars with masses from 
0.7 to 3 Mq when the energy from DM annihilation compensates 
the gravitational energy during the collapse. These equilibrium 
positions, where stars will remain for an indefinite time, are plot- 
ted for different dark matter halo densities, indicated in units of 
log(p^/GeV cm~^) at the side of each line. The grey lines are 
the classical evolutionary paths, which these stars follow before 
stopping. The DM particles are as described in Figure [T] 




3.9 3.85 3.8 3.75 3.7 3.65 3.6 3.55 



log ( T^ff I K ) 

Fig. 5. — HR diagram of a star of 1.5 M© evolving in a halo of 
DM with = 10^^ GeV cm~^. The stationary states are new 
equilibrium states where the star is powered only by energy from 
DM annihilation. These states are approximately equivalent either 
if they are reached from the collapse of the protostar (red line), or 
considering that the star evolves in a halo of DM from the ZAMS 
(green line). The grey line shows the normal HR track of a 1.5 M© 
star without DM. The DM particles are as described in Figure^ 

= 3 • 10^^ GeV cm~^ completely stops its collapse at 
the age of 50 Myr and remains forever with a stellar ra- 
dius = 1.75 R0, effective temperature Tg// = 4555 K 
and luminosity = 1.2 (Cf. Figure |4]). 

The addition of more WIMPs dramatically increases 
the central luminosity of these stars, requiring the con- 
vective core to grow in order to transport the additional 
energy to the surface layers. In this scenario, the surface 
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Fig. 6. — Tracks on the HR diagram of the collapse of stars of 
different masses until the ZAMS (grey lines), together with the 
isochrones of different ages, for DM halo densities = and 
= 3 • 10"'^'-' GeV cm~^. The DM particles are as described in 
Figure [T] 

convection zone merges with the inner core and the star 
becomes fully convective. 

Our approach on the evolution of low- mass stars within 
DM halos consists in considering th e influence of DM 
capture since the collapse of the star. 'locco et all (j2QQ8) 
did the same to study the evolution of flrst stars; they 
also found that the collapse of these stars, the so-called 
dark stars, may be stopped at an early stage if there is 
enough DM on their interior. On the other hand, an- 
other approach considered in the literature is to evolve, 
from the ZAMS, stars that were already formed with- 
out DM. T his scenario was flrst analytically estimated by 
iSalati fc'S ilk (1989) for main se quenc e stars and recently 
numerically simulated by Scot t et al.l (|2009) for low-mass 
stars. To compare the two different scenarios, we evolved 
stars using both approaches and found that they lead 
to equivalent flnal equilibrium states, even though the 
tracks followed by these stars are completely different. 
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Fig. 7. — Tracks on the HR diagram of the evolution through 
the MS (from Sgrav < IYoSt to Xc < 0.001) of stars of different 
masses (grey hues), together with the isochrones of different ages, 
for DM halo densities Px = and Px = ^ • 10^° GeV cm-^. The 
DM particles are as described in Figure [T] 

When evolved from the ZAMS, stars go back through 
the pre-main sequence phase, where they reach the same 
equilibrium states (fuelled only by the energy from DM 
annihilation) than those obtained when the collapse is 
frozen (see Figure [5j). 

4.2. Isochrones for low-mass stars evolving in DM halos 

As a synthesis of the new stellar evolution scenarios 
presented here, we show in Figures [6] and [7] the paths 
on the HR diagram followed by stars from 0.7 to 3 M© 
that form and evolve in halos of different DM densities 
p^, as well as the isochrones of different stages during 
their evolution. In Figure [6] we plotted the collapse un- 
til the ZAMS of stars that form in halos with densities 
Px = and Px = ^ - 10^° GeV cm"^. Note that, in the 
latter case, stars with masses < 1.8 Mq completely stop 
their collapse before reaching the ZAMS (as can be seen 
from the 10 Gyr isochrone). More massive stars are less 
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Fig. 8. — Effective temperature (a) and luminosity (b) of stars 
of 0.7, 1 and 1.5 as a function of the dark matter halo density, 
considering different WIMP-proton spin-dependent scattering cross 
sections. At the top horizontal axe of each figure we show at which 
distance to the galactic center are these DM densities expec ted to 
be found, following the profile of Berto ne &: MerrittI (|2005l ). The 
DM particles are as described in Figure [T] 

affected by DM; their classical evolutionary path is only 
slightly delayed in a halo with the same DM density. This 
can also be seen by looking at the 1000 Myr isochrone 
in Figure [71 In this Figure we plotted the paths of these 
stars through the MS, that is, from Sgrav < 1% St to 
Xc < 10-3. 



5. 



STELLAR DIAGNOSTIC ON THE NATURE OF DARK 
MATTER PARTICLES 



We have computed luminosity and effective tempera- 
ture of stars in all scenarios for different values of spin- 
dependent cross section (Jx,sd and DM density p^, which 
we varied from 10^ GeV cm^^ up to 10^^ GeV cm^^. 
The results are displayed in Figure [HI which also shows 
the possible location of these stars towards the centre 
of our Galaxy, where the highest densities are expected 



(|Navarro et al.l 119961 : [Salvador- Sole et al.ll2007[ ). In that 
region, the DM density may be enha nced due to the 
prese nce of a supermassive black hole (jGondolo fc Silkl 
[19991 ). We computed the DM distribution around the 
center of the Milky Way using the adiabatic contracted 
profile of Bertone & Merritt (2005). The predicted Te// 
and in Figure [8] may be used as an alternative method 
to constrain the WIMP-proton SD scattering cross sec- 
tion cr^^sD^ to help in the validation or rejection of DM 
particles' models, as well as to infer indirectly the DM 
density in the place where the star is observed. It is 
worth noting that, at the present moment, these results 
offer a qualitative picture more than an exact approach, 
due to the uncertainties in our knowledge of the inner 
region of our galaxy. In addition to the density profile 
Pxi^Gc)^ both the velocities of the star and of the DM 
particles also play an important rol e wh en studying the 
stars at the galactic center. Scott et al.[ (|2Q09 ) did pre- 
cise simulations of possible orbits of low-mass stars in 
that region, and found that only those stars with ellipti- 
cal orbits are efficient at capturing DM particles. 

In Figure [8] are plotted the Tg/ / and at such an age 
that all the stars are already in energy equilibrium. Stars 
that evolve on DM halos of low densities (Weak scenario) 
are in equilibrium in the beggining of the MS. As p^ in- 
creases, the curves mimic a slower evolution through a 
classical evolution track on the HR diagram (see Fig- 
ures [6] and [7]). For high p^ (Intermediate and Strong 
scenarios), stars are in equilibirum, powered only by the 
energy from DM annihilation. The higher the value of 
p^, the sooner the star will freeze its position on the 
HR diagram at a lower effective temperature Tg// and 
a higher luminosity (see Figure [f]). In the case of 
1 Mq star, the Tg// decreases more than 10^ K and the 
raises up to three times higher than in the MS. The 
rapid drop in temperature is related to the fact that the 
star becomes fully convective. Our result s are more con- 
servat ive than the similar ones found by [Fair bairn et al.[ 
(j2008[ ). although they predicted the same behaviour of 
Teff at lower DM densities. Probably, our underestima- 
tion could be overcome by increasing the resolution of 
our code in the very center of the star, where the energy 
from DM annihilation is produced. 

6. CONCLUSIONS 

One of the consequences of the formation of structures 
in the Universe is the creation of localised regions with 
high concentrations of DM. The formation of stars in 
such peculiar neighbourhoods should be quite different 
from the usual picture of the formation of young stars 
by gravitational collapse. A striking observational case 
of such high density DM regions is the young stellar for- 
mation regions near the super-massive black holes, lo- 
cated in the centr e of galaxies, such a s our own Milky 
Way (|Genzel et al.[ [20031 : [Krabbe et al.[[1995h . In the at- 
tempt to grasp the formation of young stars in such un- 
expected neighbourhoods, in this work we have shown 
some numerical simulations of stars with masses from 
0.7 to 3 Mq evolving within halos with high density of 
annihilating DM particles. 

We have found that the evolution of a young star can 
be affected slightly, moderately or strongly depending on 
the DM density of the host halo. Conveniently, we chose 
to classify the formation and evolution of low-mass stars 
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Fig. 9. — Classification of the different cases of stellar evolution 
within dark matter halos for stars of 0.7, 1, 1.5 and 3 Mq. Labels 
W C, I C and aS* C indicate, respectively, that these stars are 
expected to evolve in the Weak, Intermediate or Strong scenarios, 
considering halos with different DM densities and WIMP-proton 
spin-dependent scattering cross sections. The DM particles are as 
described in Figure [T] 

in three major possible evolution scenarios: Weak, In- 
termediate and Strong, which are directly related to the 
amount of DM density. The evolution of the star also de- 
pends on the scattering cross section of the DM particles. 
Figure [9] shows in which of these three stellar evolution 
scenarios should stars of 0.7, 1, 1.5 and 3 Mq evolve, for 
different values of the SD scattering cross section and the 
DM density in the halo. 

We have shown that low-mass stars in the Strong and 
Intermediate cases evolve quite differently from the clas- 
sical path of stars on the HR diagram. During their 



pre- main sequence phase, these stars will stop their col- 
lapse before reaching enough central temperature to start 
thermonuclear reactions, and will remain indefinitely in 
the same position on the HR, fuelled only by the energy 
from WIMP annihilation. In particular, stars immersed 
in high density DM halos have their effective tempera- 
ture and luminosity strongly affected due to the change 
in the energy transport in their interior (Cf. Figure [8]). 

The new data obtained by means of the near-IR in- 
strumentation allowed the ob servati on of stars in the in- 
ner parse c of our own g alaxy (jLu et al. 2009; Ghe z et aTl 
120051 : lEisenhauer et al.l ii2005). The observations have re- 
vealed a population of apparently young stars in this re- 
gion, whose current conditions seem to be unsuitable for 
star formation. This stellar population is usually consid- 
ered to be a population of old stars that have followed 
quite an atypical evolution path. If some of these stars 
were found to be low-mass stars, they would become 
candidates for this new population of stars evolving; in 
DM h alos, as initially suggested by iMoskalenko fc Wail 
([2006h . 

If found, such stars would be interesting probes of DM 
particles near super-massive black holes. Their lumi- 
nosity, or rather their excess of luminosity, attributed 
to WIMPs' burning, can be used to derive the WIMPs' 
matter density at their location. On the other hand, the 
lack of such unusual stars may provide constraints on 
WIMPs' density, WIMP-nuclei scattering and pair anni- 
hilation cross section. 

Finally, it is worth mentioning that more detailed 
studies should be done aimed at testing the validity of 
our model against new stellar observations. This paper 
should set the foundations for further works. 
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